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Abstract 

We present an algorithmic approach to calculate the quantum-noise spectral density of photocur- 
rents generated by optical fields with arbitrary discrete classical spectrum in coherent or squeezed 
states. The measurement scheme may include an arbitrary number of demodulations of the pho- 
tocurrent. Thereby, our method is applicable to the general heterodyne detection scheme which 
is implemented in many experiments. For some of these experiments, e.g. in laser-interferometric 
gravitational-wave detectors, a reliable prediction of the quantum noise of fields in coherent and 
squeezed states plays a decisive role in the design phase and detector characterization. Still, our in- 
vestigation is limited in two ways. First, we only consider coherent and squeezed states of the field 
and second, we demand that the photocurrent depends linearly on the field's vacuum amplitudes 
which means that at least one of the classical components is comparatively strong. 
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I. INTRODUCTION 



Investigations of the quantum-noise spectrum usually elaborate on the properties and 
correlations of the quantum vacuum while assuming simple classical components of the field 
. In another publication the authors investigate the quantum-noise contribution to a 
slightly more complex classical spectrum, i.e. the outcome of a heterodyne power measure- 
ment including one subsequent demodulation of the photocurrent |4( . However, the classical 
spectrum of fields inside real instruments is usually more complex, comprising many pairs of 
heterodyning sidebands 5|, |6j. In addition, one may be interested in measurement schemes 
which include more than one demodulation of the photocurrent. In this paper, we extend 
previous analyses by allowing for an arbitrary number of discrete classical components and 
an arbitrary number of demodulations. Furthermore, our approach is algorithmic which 
means that the results can be implemented straightforwardly in the code of simulators of 
quantum-noise spectra. 

In Sec. HIl we introduce our notational conventions and calculate spectral densities for 
power measurements of fields in coherent states. The complexity is gradually increased, 
starting with a measurement of discrete components without demodulation and ending with 
multi-component fields including arbitrarily many demodulations. These results are further 
generalized in Sec. [ITT] to include fields in squeezed states. For this purpose, the representa- 
tion of the photocurrent is slightly but essentially modified to take account of the intricate 
sideband correlations occurring in squeezed fields. Finally in Sec. HV] the algorithm is applied 
to investigate the photocurrent noise spectrum for a specific multi-component, squeezed field 
configuration which is meant to clarify the abstract approach outlined in previous sections. 



II. QUANTUM-NOISE SPECTRUM OF COHERENT FIELDS 

Let us first calculate the quantum-noise spectral density of a coherent field which is 
determined by a single classical amplitude c at frequency / . Denoting the quantum vacuum 
noise amplitudes by q(f), the electric field can be written 

E(t) = c e~ 27Tifot + c*e 27Tifot 

OO 

+ Jdf (q(f)^ ift + 9 f (/)e 27ri/t ) ■ 
o 
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Factoring out the oscillating phase, exp(— 2ni f t), of the classical field, the quantum vacuum 
integration is carried out over sideband frequencies F = f — fo- The range of sideband 
frequencies is restricted since we consider a limited measurement bandwidth: F G [-B, B] 
(B <C fo). We also split the field into two Hermitian conjugate parts E^ + \t), E^(t) and 
so the positive-frequency field is 

E ( +\t) = e" 27ri/o ' I c + J dFq(f + F)e~ 2 * iFt j . (2) 



-B 



When talking of spectral densities, we mean the power spectral density of the photocurrent 
after demodulations. Ideally, the photocurrent I(t) is proportional to the power of the field. 
For real photodiodes, this does not have to be true for arbitrarily large detection bandwidths. 
In the simplest case the current spectral density has to be multiplied with a factor which 
accounts for the frequency dependent response of the photodiode. Let us assume that B is 
small enough so that I oc P for all measured frequencies. The power of the field Eq. (j2J) 
averaged over a period r <C c/B is given by 



P(t) = M-)(t)M+\t) 



B 



l -\c»\ 2 + J dFc*q(fo + F)e- 2 ^ Ft + h.c. 



(3) 
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Contributions quadratic in the vacuum amplitudes are neglected which is a valid approxi- 
mation whenever Pq = \cq\ 2 ^> hfo/r , where r denotes the measurement time. It follows 
that the optical quantum noise of the photocurrent 7^ M (t) is proportional to 

B 

J QM (t) oc J dFc*q(f + F)e~ 2 ™ Ft + h.c. (4) 

-B 

The proportionality factor \ between current and light power — the photodiode 
responsivity — typically assumes values \ = 0.2 — 0.7 A/W at / ~ 3 • 10 14 Hz 7]. Since 
subsequent calculations are based on a frequency independent responsivity, there is no real 
distinction between light-power and photocurrent power spectral densities except for a factor 
of x 2 . Throughout this paper the explicit dependence of the photocurrent on x is omitted. 

Vacuum amplitudes like any classical amplitudes of stationary noise are numerically ill- 
defined (the Fourier transform of stationary noise being infinite) and serve exclusively as 
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algebraically meaningful quantities. Nevertheless, the power spectral density of stationary 
noise is well-defined and there exists a simple relationship between stationary vacuum noise 
amplitudes q and its (single-sided) noise spectral densities S q 

2 \ (5) 

= 2*(/ - /') h f- 

The corresponding quantum noise spectral density of the photocurrent Eq. (J4]) which results 
from the power measurement of a coherent field reads 

S? M (F) = \c \ 2 -S q (f ) = P -hf . (6) 

The invalidity of the latter equation for squeezed states of the field and the fact that the 
spectral density does not depend on the sideband frequency F is clarified in Sec. IIHI (see 
Eq. (|28|) ). At this point, we only wish to draw attention to the apparently simple algorithm 
for coherent fields which leads to Eq. ([6]) starting from Eq. ([2]): square the classical amplitude 
and multiply the result by h Jo- 
in fact, the algorithm stays valid for a wider class of power measurements. Consider a 
set of classical components with amplitudes c; = c(/i), i G {0,1, . . . , N — 1}. Like before, 
we factor out the phase of the component cq which formally converts all optical oscillations 
into sideband oscillations with respect to the reference frequency /q. Defining F[ = j\ — / , 
the positive-frequency field is cast into the form 

N— 1 Fi+B 

EW (t) = e~ 27Tif0t + J dFq(f + F)e~ 2 * iFt J . (7) 

i=0 Fi-B 

We demand that P\ = |c;| 2 ^> hfo/r is valid for all classical components. If that condition did 
not hold for some classical components (typically for signal sidebands), then the respective 
amplitudes would not influence the quantum-noise power spectrum as long as at least one 
high power component exists and so for our purpose we may safely neglect weak components. 
A power measurement yields a photocurrent whose quantum noise is determined by 

N-l Fi ~t B 



fQM( t) = J2 I dFc*q(f + F)e- 2 ^ F - F ^ + h.c. 

(8) 

N-l ° 

= J2 <iFc*q(f + F + F^ Ft + h.c. 
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Provided that all intervals [Fi — B,F{ + B] are mutually disjoint, the quantum-noise power 
spectral density of the photocurrent is given by 

N-1 N-1 

S? M (F) = |Q| 2 • S q (/o + Ft) = h ■ (fo + Fi) ■ ^ (9) 

i=0 i=0 

Unfortunately, there are no more simple example cases which could be presented. Let us 
now skip to the next section and introduce demodulations of photocurrents. 



A. A Single Demodulation 

There are two well-known demodulation techniques. Either the photocurrent is multiplied 
by a harmonic function I(t) ■ cos(27r.D -t+<f>) or one rectifies the photocurrent which means 
the final output is \I(t)\. In practice, the low-frequency spectra [B <C Fi) drawn from 
these two outputs differ by a constant factor, the spectrum of the harmonic demodulation 
being smaller by a factor of (7r/4) 2 . Anyway, the spectrum of the harmonically demodulated 
current is calculable by much simpler algebra. In addition, we want to study measurement 
schemes which implement multiple demodulations characterized by a set of demodulation 
frequencies D x and demodulation phases 0;. For these two reasons, we consider harmonic 
demodulations throughout this paper. Let us perform a single demodulation of the current 
Eq. © 

N-i - D + B 

N-1 D + B ( 10 ) 



+ S &Fclq{h + F + F t )e-^ F - D > 



1=0 D-B 

+ h.c. 

or shifting the frequency range, the current noise finally assumes the form 

N-1 B 

I QM (t) = ^Yl I dFe-^ Ft [c:e-^q(fo + F + F t -D) + c*e^q(f + F + F l+ (11) 



i=0 



-B 



Now, if all frequency intervals [Fi — D — B, Fi — D + B] , [Fi + D — B, Fi + D + B] are mutually 
disjoint, then one obtains exactly the same spectral density as in Eq. j9}. However, this 
situation is uncommon in real experiments. 



It is time to introduce a graphical auxiliary in order to understand what is happening. 
For simplicity we start with two classical components — -a carrier and a subcarrier — at fre- 
quencies F = Hz and F% = 30 MHz. The photocurrent is demodulated with D = 15 MHz 
and the detection bandwidth is limited to B — 1000 Hz. This choice of frequencies also 
corresponds to a very common situation in optical experiments, namely the heterodyne 
measurement where the two classical components are generated by a 15 MHz modulation of 
a carrier field which oscillates at high optical frequencies (fo ~ 10 15 Hz). We point out that 
B -C \Fq ± D\, \Fx ± D\, \Fq ± F\\ which significantly simplifies the problem. The four fre- 
quency values F o = F Q — D, F i = F Q + D, F w = F\ — D and F\\ = F\ + D are marked on a 
frequency axis (see Fig. [TJ) and collected inside a matrix F(N = 2, 2) = {F 00 , F 01 ; F 10 , F u }. 
Since two intervals identically overlap, each vacuum amplitude at frequencies inside that 



Fqo Fqi — Fio Fn 



FIG. 1: Representation of the relevant frequency ranges for a classical heterodyne power measure- 
ment. The black boxes indicate the detection bandwidth. 

interval enters twice into the spectral density calculation. These two contributions have to 
be added coherently before taking the absolute square. The expression Eq. (flT]) for the 
current tells us that all vacuum amplitudes q(fo + F+ 15 MHz) have to be added coherently. 
In conclusion, the spectral density in this particular example turns out to be 

S? M {F) = ~ (h ■ (f - 15 MHz) • P + h ■ (f + 15 MHz) • lege" 1 * + c*e^| 2 

+ M/ + 45MHz)-P 1 ) (12) 
^^(P + P 1 + |c*e-^ + C ^| 2 ). 

Setting ci = 0, one may wonder where half the noise power has gone when comparing 
the result with Eq. Basically, the answer is that by demodulating we multiply the 

photocurrent with a function which, in terms of power, has a gain of 1/2. 

Henceforth, we will always assume that overlapping frequency ranges represented by black 
boxes in Fig. [1] overlap completely, but never partially. The overlap condition guarantees 
that whenever two frequencies inside the frequency matrix T do not coincide, then the two 



6 



respective detection ranges do not overlap and one does not have to worry about coherent 
summation of the respective amplitudes. This is certainly a reasonable demand for a first 
approach to an algorithmic realization of the calculation. 

A discussion of the result Eq. f[T2"j) focussing on signal-to-noise ratios can be found in 
4j. We are now prepared to calculate the spectral density for measurements with a single 
demodulation and an arbitrary number of classical components. It is not possible to give 
explicit results, because as we have seen these depend on the chosen set of field and demod- 
ulation frequencies. Our focus lies on extending the algorithm which leads to the spectral 
density. Consider N classical components c\ at frequencies F x and a single demodulation of 
the photocurrent. The first step is to calculate the frequency matrix 



( F F \ 

-^00 ^01 



(13) 



where F\§ = F\ — D and Fn = F\ + D. The second step is to find coinciding frequencies of 
the matrix. Use the matrix indices of these pairs, e.g. (ni,di) and (712,^2), to calculate the 
contribution to the spectral density as follows 

S? M (F, K dO, (n 2 , d 2 )) = h -^ + Fnidl) ■ K 1 e (_1) * +1 * + 4 2 e ( " 1)d2+li 1 2 (14) 

All remaining unique frequencies F n d contribute according to 

Sr(F,(n,d)) = h - {fo + Fnd) \c n f. (15) 

In most experiments /q ^> F n( ^ and the quantum vacuum energies can be approximated by 
hfo- Finally, one has to sum up all these contributions. 

Let us summarize our preliminary results as a list to be processed when calculating the 
current noise spectral density: 



1. Calculate the matrix T according to Eq. ([13]) . 

2. Collect index pairs (n, d) of coinciding frequencies F n d inside T . 

3. Collect contributions to the current noise spectral density from unique frequencies 
which are determined by Eq. ffT5|) . 



4. Collect contributions to the current noise spectral density from coinciding frequencies 
which are determined by Eq. (jl4j) . 

5. Sum up all contributions. 

This algorithm needs an extension in order to account for multiple demodulations. This 
task is accomplished in the next section. 

B. M Demodulations 

Methodically, increasing the number of demodulations means to increase the number 
or range of indices in the previous calculations. The first step of the algorithm is always 
to collect all relevant frequencies inside a matrix. Again, we consider a set of iV classical 
components c\ at frequencies F\. The demodulations are determined by M demodulation 
frequencies and phases 4>- x . The final output of the measurement is given by 

M-1 2 M /M-1 \ 

I(t) ■ J] «*(2*A + = lit) • E cos E (-1) Li/2M " J_1J (2vr Z)j • t + 0,) (16) 

i=0 i=l \ j=0 / 

with \_x\ denoting the decimal truncating floor function. Evaluating Eq. ffTBT) for a few small 
number of demodulations, one recognizes that the arguments of the harmonic functions 
consist of all possible ±-combinations of modulation frequencies and phases. Since the 
complete sum of these harmonic functions is multiplied with the photocurrent, it should 
be obvious that each row i of the frequency matrix T contains all frequencies of the form 
Fi ± D ± . . . ± Dm-\- In which order should these 2 M frequencies appear? One may argue 
that the order does not matter in principle. However, a specific combination of frequencies 
is associated with an analogous combination of demodulation phases which then appears 
in formulas like Eq. ffT4l) . If two frequencies of T turn out to be equal, then the two 
pairs of matrix indices have to encode the ±-combination which determines the respective 
combination of phases. Algorithmic tractability of the problem requires a good sorting 
scheme of these frequencies. We propose a sorting scheme which is derived from a tree 
structure like that in Fig. [2J The frequencies themselves do not obey F i0 < Fn < . . . since we 
do not assume a magnitude-sorted vector of demodulation frequencies. Now, the remaining 
problem is to derive from each column index the respective combination of demodulation 
phases. If each row vector of T is tree sorted then one first converts the decimal column 
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+Dq 
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-D\ +D 1 -D x +D X 

-D 2 -D 2 +D 2 -D 2 +D 2 -A 
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F\q Fix Fj2 Fr Fj4 F;5 F;6 -F; 
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FIG. 2: Sorting scheme for three demodulation frequencies. 



index j into a binary number (j)un = [bobib 2 bz . . . and then calculates a total phase 

according to 

M-l 

M = -Y.(- 1 ^^ ( 17 ) 



8=0 



This phase is to be used in order to calculate contributions to the quantum noise spectral 
density when C frequencies of T coincide. If the respective pairs of indices are (n\,dj), 
j £ {0, .., C — 1}, then the spectral density reads 



h ■ (Jo + F nodo ) 
4M 



c-i 

E 

3=0 



(18) 



The contribution from unique frequencies F n< j is given by 

SrWn,d)) = h - {f l+ Fnd) \cn\ 2 . 



(19) 



Finally, let us apply the algorithm to the case with N = 2 classical components and M = 3 
demodulations. The respective frequency matrix is written 



^(2,8) 



(20) 



FqQ Fqi Fq 2 F03 F04 F05 Fqq Fqj 

F10 Fu F\2 Fx3 F 14 F 15 F 16 F 17 
Coincidences are found say between frequencies F 05 and Fu. The binary representation of 
the first column index is (5)bi n = 101 and (l)bin = 001 for the second. According to Eq. i fPTt) . 
one obtains 0(5) = </> — 4>i + 02 and 0(1) = — <fi — <f>i + 4> 2 . The respective spectral density 
is given by 

~ n ™-- ' - h-(fn + Fu) ,, 

Gge -|- c^e 



S? M (F, (0, 5), (1, 1)) = M/o + fli) , c * e i0(5) + c * e ^(i),2 (21) 



64 

Contributions coming from unique frequencies are calculated as usual. The first part of 
this paper is finished. We have treated the calculation of quantum noise when the field is 



coherent which entails that the noise is time stationary. In fact, the algorithm is valid for 
any type of time-stationary noise including technical laser noise when phase and amplitude 
noise are not correlated. In that case, one has to substitute the quantum vacuum energies 
in Eqs. (I18 p &; (|19]) by another noise spectral density which characterizes the power of the 
technical noise at the photodiode. The next step is to take into account the intricate 
correlation between quantum noise amplitudes at different frequencies due to squeezing 
or in other words, due to amplitude-phase correlations. 



III. QUANTUM-NOISE SPECTRUM OF SQUEEZED FIELDS 

A widely applied mechanism which correlates vacuum noise amplitudes at different fre- 
quencies is the generation of squeezed fields. Squeezed fields are formed in nonlinear crystals 
Q and theory predicts that pondero.not.ve Sq nee,n g occnrs when light is reflected from 
suspended mirrors [10]. In order to understand the meaning of squeezing, one has to know 
that correlations are built up between sideband frequencies F with respect to a reference 
frequency f\. The squeezing transformation of fields is characterized by a squeezing factor r 
which quantifies the strength of correlations between different frequencies and a squeezing 
phase 0. Perfect correlation corresponds to a squeezing factor r = oo. In practice, the 
reference frequency is realized by means of a single classical component which is named seed 
or carrier field depending on whether the squeezing is generated by crystals or ponderomo- 
tively. Now, if = one can show that the correlation between sidebands diminishes the 
quantum amplitude noise of the classical field and for <p = tt/2, the quantum phase noise 
is decreased. Interpretation of the squeezing phase for intermediate values requires a more 
sophisticated representation of fields jg|. In principle, the reference frequency does not have 
to be related to a classical field. In that case, the squeezing phase has no ad hoc interpre- 
tation. Denoting amplitudes of squeezed fields by s and amplitudes of coherent fields by q, 
the squeezing transformation is governed by 8| 

s(fi ±F) = cosh(r) • q(fi ± F) + sinh(r)e 2i< ^ • q\f { T F). (22) 

How does this squeezing transformation relate to our previous investigation of quantum noise 
spectral densities? To find an answer, we have to return to Eq. (jSj). That equation determines 
the time-dependent photocurrent where F, the former optical sideband frequencies, now 
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become true frequencies of the current spectrum. Similarly to expansions of electric fields, 
we had better conform to a strict decomposition into positive and negative frequencies 

N-l B . 

jQM (t) = Y t j dFe- 2 * iFt (c*s(f + F + F l ) + Cl s\h -F + F<)) + h.c. (23) 

where the vacuum amplitudes have been renamed to indicate the possibility of squeezing. 
Thereby, the impact of Eq. (|22|) on the power spectral density of the photocurrent becomes 
clear. Two amplitudes at different frequencies of the optical vacuum are added to form a 
single amplitude of the photocurrent. If these two amplitudes are uncorrelated then the 
spectral density is calculated as before making implicit use of the identity 

S(qi + ga) = S(qx) + S(q 2 ), q\, q 2 uncorrelated. (24) 

The optical field may exhibit correlations due to squeezing and the latter equation can no 
longer be applied. So the remaining problem is the calculation of the modified spectral 
density depending on the squeezing factor and phase. There exist two possibilities. Either 
squeezing is generated whose reference frequency coincides with one of the frequencies fo + 
Fi which appear in the current expansion Eq. (|23|) and the squeezing factor is significant 
only within a frequency range comparable to the detection bandwidth B, or the squeezing 
violates either of these two conditions. We start with an investigation of the case when 
both conditions are fulfilled which is in some sense the only expedient one. It is then also 
straightforward to treat fields which are squeezed at some or all of the frequencies fo+Fj with 
potentially different factors and phases. The correction is derived from the quantum-noise 
amplitude of the photocurrent; 

p(F, c*, F ; ) = c*s(f + F + F i ) + c^(f -F + FJ. (25) 

Before or after squeezing, the field may be subject to linear (frequency preserving) transfor- 
mations like propagations or reflections from fixed mirrors. If some mirrors are suspended 
then already squeezed fields experience further squeezing. Consequently, a more generic case 
is considered here described by the following pair of transformations: 

s(f i + F)=t ao (f i + F)q{f i + F)+t m .(f i -F)q1(f l -F) 

(26) 

— F) — t u (/, + F)q(fi + F) + t n (fi - F)q\}\ - F). 
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The two transfer functions tio, toi m &P input amplitudes to output amplitudes at the pho- 
todiode with mutually different frequencies, whereas too, t\\ comprise nonlinear and linear 
transfers as mentioned above. The input vacuum is coherent and therefore all amplitudes q 
of the input vacuum are uncorrelated. Inserting the last equation into the amplitude of the 
photocurrent and defining tj nn (±) = t mn (/o ± F + F{), one obtains 

p(F, cf, Ft) = [c*4(+) + c4 (+)] • q{f + F + F t ) 

+ [c:4(-) + Qt 1 11 (-)]-gt( /o -F + F l ). 

The power spectral density associated with the noise amplitude p(F, c*, F{) is given by 



S(p(F,clF i )) = h -( f ° + Fi) 



|c*4(+) + ^i (+)l 2 - (i 



fo + Fi 

+ic*4H+c4i(-)r (i- 



(28) 

F 



fo + F 



First, this equation explains the observation that the spectral density of photocurrents when 
measuring coherent fields is white (i.e. frequency independent). For coherent fields (£qi = 
t w = 0, |too| 2 = |^n| 2 — 1)) the two absolute squares can be substituted by |cj| 2 and the 
dependence on F cancels. Since in most practical situations /o ^> F, it is also common to 
omit the frequency dependence of the spectrum for squeezed states. 

The next step is to include demodulations of the photocurrent. In terms of the amplitudes 
defined in Eq. ( |271) . a singly demodulated photocurrent, Eq. ( fTTT) . assumes the form 

/Q M (t) = ~J2 ^Fe- 2 ^ Fd (p(F, c*e-^, F { - D) +p(F, c^, F + D)J + h.c. (29) 







Hereupon, a simple analogy argument leads to the full-fledged multiple demodulation power 
spectral density of the photocurrent including squeezed states of the field. If C frequencies 
with indices (nj, c?j), j G {0, .., C — 1} of the frequency matrix T coincide, then the respective 
contribution to the power spectral density is determined by (compare with Eq. ([TBI) ) 
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sFV.^.d,)}) = ^-s y£p^ci^\F nodo ) 

\j=0 



M/o + F, 



no do i 



2 • 4 M 



c-i 



+ 



3=0 



C-1 



' : /,H,(./u • 

F 



F nodo ) + c^e-^hoifo + F + F nodo ) 



1 + 



fo + F 



n d Q 



(</ m) toi(fo -F + F nod0 ) + c^e-^tnC/o - F + F„ odo ; 
i=o 



/o + F 



no do 



and unique frequencies F nd contribute with 



S? M (F,(n,d)) = ^ ■ S (p(F,c n ,F nd )) 



(30) 



(31) 



The spectral density on the right-hand side is determined by Eq. (l28l) . 

We conclude this section with a brief discussion of "off-centered" squeezing, i.e. some of 
the squeezing reference frequencies of the field do not coincide with any of the frequencies F nc i 
or that squeezing factors r decay over frequency ranges which are comparable to F nc j 3> B. 
We leave a detailed investigation of both problems for the future. One reason is that in each 
case amplitudes of the measured field at the photodiode may depend on input amplitudes at 
more than two frequencies. In other words, one squeezing process may correlate amplitudes 
at say j\ = 10 15 Hz + 50 MHz and f% = 10 15 Hz + 30 MHz, another squeezing process then 
correlates amplitudes at /a and fs = 10 15 Hz — 30 MHz. The output field at fi depends on 
input frequencies fx, f 2 and / 3 . Multiple correlations are prevented by demanding that the 
frequency difference of correlated amplitudes is sufficiently small so that different correlations 
occur at well separated parts of the spectrum. If that condition is fulfilled but "off-centered" 
squeezing is still allowed, then one has to take greater care when calculating the transfer 
functions ty. What if the squeezing reference frequency happens to be at fo + F + 5-B which 
is close enough to the measured frequency range [fo — B, /o + B\ to exhibit some influence 
on its vacuum noise amplitudes? The answer is derived from Eq. (1221) . Let us assume that 
the squeezing factor and phase are constant over [/ — B, f + B}. Then, provided that 
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B -C fo, the photocurrent spectral density, Eq. (jSJ), is modified according to 

Sj QM (F) = P • hf ■ (cosh 2 (r) + sinh 2 (r)) 
= P -/i/ -cosh(2r). 

which means the quantum noise spectral density is necessarily increased! Therefore, any 
detector whose sensitivity is limited by optical quantum noise should avoid off-centered 
squeezing. 

Again, the previous analysis is not restricted to quantum noise spectra. The results 
are valid for any classical, technical noise which is coherent or exhibits amplitude-phase 
correlations which can be described by transformations like Eq. ( |26l) . One may say that 
these transformations embody the simplest kind of amplitude-phase correlations. 



IV. EXEMPLARY CALCULATION 



^(4,2) 



(33) 



Concluding the paper with an explicit application of our results should be helpful. We 
consider the following case: four classical components, one demodulation and squeezing 
centered at one of the components F nd of the frequency matrix T. The N = 4 classical 
components have frequencies F = Hz, Fx = 100 MHz, F 2 = 130 MHz and F 3 = 230 MHz. 
These sideband frequencies are defined with respect to a large optical frequency fo ~ 10 15 Hz. 
The photocurrent is demodulated with Do = 15 MHz. These parameters determine the 
frequency matrix 

/-15 MHz 15 MHz \ 
85 MHz 115 MHz 
115 MHz 145 MHz 
\ 215 MHz 245 MHz/ 

We assume squeezing centered around frequency Fn = F20 = 115 MHz which extends 
locally over frequencies comparable to the detection bandwidth B = 1000 Hz, all other 
frequencies contribute coherent vacuum noise. We already stated that two components of 
the matrix T are equal and the respective frequency value coincides with the squeezing 
reference frequency. Let us first evaluate the contributions from all unique frequencies. The 
corresponding vacuum fields are coherent and so the spectral density Eq. (128!) simplifies to 

h{f + F nd )\c n \ 2 

Hfo + F nd )P n . 



ricoh 
D nd 



(34) 
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These spectral densities have to be inserted into Eq. (I3T]) and summed up for all unique 
frequencies 

h 



Snm(F) = - ((2/o + ^oo + F 01 )P + (fo + F 10 )P 1 
+ (fo + F 21 )P 2 + (2/o + F 30 + F 31 )P 3 



(35) 



The remaining problem is to calculate the correlated spectral density at frequency Fn 
F 2 o = 115 MHz. For simplicity we assume that the nonlinear transfer is pure squeezing 



*oo(/o + F + F u ) = cosh(r) 
hiifo-F + Fn) = sinh(r)e i2 * 
ti (fo + F + Fn) = sinh(r)e- i2 * 
tn(fo-F + F 11 ) = cosh(r). 

Next, we expand the sums in Eq. fl30l) and substitute all known parameter values 



(36) 



^ M (F,{(1,1),(2,0)}) = ^ 



+ 



c*^ + c* 2 e' i<t>0 ) cosh(r) + (cie^ + c 2 e i4>0 ) smli^e" 2 ^ 
'c\^° + c* 2 e-^) sinh(r)e 2i * + (cie -1 * + c 2 e i4>0 ) cosh(r 



(37) 

All vacuum energies are approximated by hf and 0o denotes the demodulation phase. 
There are at least two experimentally adjustable phases, the demodulation phase 0o an d the 
squeezing phase 0. What is the minimum of the spectral density depending on these two 
phases? Defining a = arg(cie~ 1</io + c 2 e^°), Eq. (1371) assumes the form 



S? U (F,{(l^),(2,0)}) = ^\c* 1 e i *° + c*e 



cosh(2r) + sinh(2r)e 



2i(4>-a) 



(38) 



Minimization with respect to the squeezing phase is trivial. Setting opt = tt/2 + a and 
further defining Aa = arg(c2) — arg(ci), the minimized spectral density contribution from 
the squeezed part of the spectrum simplifies to 



S? M (F, {(1, 1), (2, 0)}) = M( Pl + p 2 + 2^P 2 cos(20 o + Aa)) • e~ 2r . (39) 
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Finally, we add this result to Eq. (1351) which yields 



^ (Po + i(l + e" 2r ')Pi + \{l + e' 2r )P 2 - v/PP • e~ 2r + R 



(40) 



In the last step, we have minimized the noise power by inserting the optimal demodulation 
frequency 0°j Pt = (n — Aa)/2. The reader who is exclusively interested in minimized spectral 
densities may easily generalize this final result including squeezing at different frequencies 
and an arbitrary number of classical components. The optimization procedure partly relies 
on the fact that the squeezing factor and phase are frequency independent. In general, ex- 
perimental realization of the smallest possible noise spectral density requires further degrees 
of freedom which are incorporated into the transfer functions ty, i.e. the light has to be 
filtered before the power measurement Q, ll|. One should keep in mind that, per se, a min- 
imal noise spectral density does not have to be optimal in terms of detector sensitivity. The 
reason is that noise minimization simultaneously affects the measured power of the signal. 
Sensitivity optimization severely depends on the detector topology 11[ . It was shown in 4j 
that the sensitivity optimizing demodulation frequency is 0[] pt = —Aa/2 which goes along 
with a maximized noise power contribution from Eq. (1391) . 

V. CONCLUSION 

We have presented explicit formulas which govern the power spectral density of pho- 
tocurrents generated by power measurements of coherent and squeezed fields. We have 
also furnished an appropriate algorithm which can be easily implemented in quantum-noise 
simulations. The algorithm is based on a few limitations concerning the classical spectrum 
and the squeezed spectrum of the field. However these limitations are modest and 
probably insignificant for most experiments. Our results provide a long-sought-for extension 



of the Schottky formula [12|, ll3| to squeezed photon statistics with multiple classical 
components. Of special importance is that the shot noise spectrum may be calculated 
for any of the currently operating and next generation interferometric gravitational wave 
detectors. Theoretically, the method can be generalized in two ways which we consider 
to be algorithmically tractable in principle. First, the transfer functions between different 
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frequencies - which correspond to classical nonlinearities - may couple more than three 
amplitudes at different frequencies. Thereby, multiple squeezing centered around different 
frequencies with overlapping ranges of non-vanishing squeezing factor could be described. 
Second, one may want to give up the overlap condition between different ranges of detected 
field amplitudes. In that case, correlated contributions to the final current spectral density 
have to be revealed by means of a more elaborate scheme which calculates the intersection 
boundaries of a partial overlap. 
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